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Optical Feshbach resonance is capable of inducing spatially varying interactions in ultra-cold atoms. Its 
applications to pancake-shaped clouds of bosons and fermions enable one to study several fresh phenomena. 
We examine possibilities of inducing counter-intuitive structures such as creating a superfluid enclave inside a 
Mott insulator for bosons and a normal-gas core enclosed by a superfluid shell for fermions. We discuss feasible 
experimental setups and signatures of those interesting structures, which can be very different from common 
structures observed in experiments so far. While a superfluid enclave in a Mott insulator can be useful for 
constructing atomic devices for atomtronics, the properties of the superconducting islands observed in scanning- 
tunneling microscopy of heavily underdoped high-temperature superconductors may be studied with cold Fermi 
gases with spatially varying attraction. 

PACS numbers: 32.8().Qk, 74.20.Fg, 67.85.Lm, 72.80.-r 



I. INTRODUCTION 

Recent progress of experiments using ultra-cold atoms al- 
lows for simulations of complex quantum systems which can 
be formidable challenges to conventional setups lUlll. Tun- 
able atomic interactions and trapping potentials have brought 
us exciting phenomena such as BCS-Bose-Einstein conden- 
sation (BEC) crossover of fermionic superfluids and Mott 
insulator-superfluid transition of bosons |1|. More exotic 
phenomena such as heteronuclear superfluids |3| or d-wave 
superfluids |4| closely related to nuclear physics or high- 
temperature superconductors (HTSCs) have been proposed. 
Here we show that recent development of optical Feshbach 
resonance (OFR) |l5]-|8l which can control the scattering length 
in a selected region will bring another advantage of using 
ultra-cold atoms to demonstrate even more exciting quantum 
phenomena induced by controllable spatially varying interac- 
tions. Moreover, we argue that this technique is a key compo- 
nent in providing crucial links to other fields such as quantum 
simulations of electronic devices and HTSCs. 

Atomtronics [9J, whose goal is to simulate electronic de- 
vices using ultra-cold atoms, is a joint venture of several 
fields. A major challenge in atomtronics is to produce a de- 
signed pattern of atomic devices. In addition to engineered 
potentials, one can take advantage of spatially varying inter- 
actions, which can induce a superfluid region inside an insula- 
tor or vice versa as we will demonstrate here. A combination 
of the two techniques, spatially varying potential and interac- 
tions, will help realize integrated atomic circuits, the analog 
of integrate circuits of electronics. 

An important property of spatially varying interactions gen- 
erated by OFR is its capability of locally tuning the local den- 
sity for a system in equilibrium. For bosons in the ground 
state this can be seen from Gross-Pitaevskii equation or Bose- 
Hubbard model, where the local chemical potential (defined 
as the difference between the equilibrium chemical potential 
and the local confining potential) depends on the local interac- 
tion strength 1 1 1 . For fermions with attractive interactions one 
can understand this from BCS-Leggett theory of BCS-BEC 
crossover of Fermi gases |10|, where the fermion chemical 



potential ^ / decreases monotonically as the attractive interac- 
tion increases. 

For a trapped atomic gas, the local density approximation 
(LDA) suggests that the chemical potential has to balance the 
trap potential. In the presence of spatially varying interac- 
tions, the local chemical potential includes effects from both 
trap potential and the controlled interactions. As a conse- 
quence, the density profiles can change dramatically when the 
strength of interactions changes from one region to another. 

A major goal of this paper is to show that counter-intuitive 
structures can arise in well studied systems such as Bose- 
Hubbard model and BCS-BEC crossover. For example, a 
normal Fermi gas can coexist with a Fermi superfluid, but 
they have to arrange themselves to balance /i/ so that there 
is no net mass current. For a Fermi gas with uniform attrac- 
tive interactions in a harmonic trap, the superfluid phase re- 
sides in the center of the trap, where the density is higher, 
to take advantage of its condensation energy. In contrast, by 
suitably shielding the attractive interactions at the trap center, 
it is possible to have a ground state with a normal-gas core 
encircled by a superfluid shell. We note that for equal-mass 
fermions with uniform interactions, the superfluid always re- 
sides at the trap center even in the presence of population im- 
balance l(3| nn . An inverted structure has been proposed in 
the ground state of heteronuclear Fermi gases with strongly 
attractive interactions |3 1 but has not been observed in exper- 
iments yet. Here we argue that a normal-gas core can exist 
in the ground state without introducing two different masses 
of fermions with mass imbalance if one can spatially vary the 
interactions using OFR. 

Another important motivation for implementing tunable 
spatially varying interactions comes from the fact that many 
interesting condensed matter systems are intrinsically inho- 
mogeneous and the effective interactions thus change spa- 
tially. For example, superconducting islands immersed in 
a non-superconducting background have been reported in 
scanning-tunneling microscopy images of heavily underdoped 
cuprate HTSCs |12|. To help understand the mechanism of 
HTSCs, experiments on ultra-cold atoms should include such 
inhomogeneity effects and study how significantly those ef- 
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fects suppress the transition temperature. 

There are other techniques for manipulating locaHzed re- 
gions in an atomic cloud. For example, a spatial light modula- 
tor (SLM) 1 13 1 can generate time-dependent patterns of trap- 
ping potentials for atoms so it is usually used as optical tweez- 
ers for atoms. There are experiments demonstrating splitting 
of BEC [14] or moving an atomic cloud |15|. Thus SLMs 
are more suitable for studying dynamic properties. In this pa- 
per we focus on static configurations that can be generated by 
OFR, but together with the SLM one may further study dy- 
namic phenomena. 

Before presenting our theoretical studies, we briefly sum- 
marize relevant experimental studies on OFR. In bosonic BEC 
of ^^"^Yb it has been demonstrated experimentally |8| that a 
modulation of the scattering length with alternating regions of 
width ^ 278nm can be generated by OFR. Such technique can 
be applied to fermionic ^^^Yb and ^^^Yb as well, which have 
been cooled down to 0.46Tf and 0.54rp [16|. Those temper- 
atures are not far from the superfluid transition temperature of 
unitary Fermi gases, which is estimated to be O.lSTTp ifTTl . 
Moreover, one can introduce bosonic ^^^Yb or ^^"^Yb into the 
fermionic system and the boson-fermion or boson-boson in- 
teractions can all be tuned by OFR ifTSl . Adding to the ex- 
citement is the successful realization of OFR in bosonic ^'^Sr 
reported in Ref. [19|. Thus our predictions not only provide 
future directions but also are readily testable as experiments 
progress. 



II. SUPERFLUID ENCLAVES IN A MOTT INSULATOR 

We first show how spatially varying interactions can induce 
a superfluid enclave inside an insulating region. Here we con- 
sider Bose-Hubbard model (see Ref. 1 1 1 for a review) describ- 
ing bosons on optical lattices in quasi two dimension (2D). 
The Hamiltonian is 

Hbh = -t'^blbj - fi'^ni + ^'^ni{ni - 1). (1) 

(ij) * * 

Here bi (bj) is the boson annihilation (creation) operator, (ij) 
denotes nearest neighbors, t and U are the hopping coefficient 
and onsite repulsive coupling constant. There is also a back- 
ground harmonic trap potential Vtrap = where 
LJtT is the trap frequency. The Mott insulator (MI) - superfluid 
transition in quasi 2D has been demonstrated experimentally 
II20I . When the onsite repulsion U is strong (weak), bosons 
localize (delocalize) and form a Mott insulator (superfluid). 
Figure [TJ a) shows the phase diagram of the n — 1 MI and the 
superfluid phase close to it (see Ref. |[T] for details), where n 
is the filling factor 

In the presence of a harmonic trap, the chemical potential 
decreases from the trap center to the edge as shown by the 
path labeled (b) on the phase diagram. Thus the n — 1 Ml 
resides at the trap center enclosed by a superfluid at the edge. 
The structure is schematically shown in Fig. [TJb). Now we 
consider tuning U using OFR because U is controlled by the 
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Figure 1. The phase diagram of 2D Ml-superfluid transition (a). 
The two paths (b) and (c) show two possible phase structures from 
the trap center to the edge in a harmonic trap following LDA with 
fi{r) — ^{r = 0) — Vtrapir). For uniform interactions, a MI core 
is enclosed by a superfluid (SF) as shown in (b), which can be un- 
derstood by following the path "(b)". When a shield blocks the en- 
hancement of repulsion induced by OFR, the shielded region remains 
a SF, as shown in (c). Due to a sudden change in U and continuity 
in /i, there are jumps in both directions in (c) 1211 . A proposed ex- 
perimental setup is shown in (d), where we emphasize that the shield 
needs not be above the trap center. The structure of (c) can be un- 
derstood by following the path "(c)". Here t and U are the hopping 
coefficient and onsite repulsion. Possible diffraction patterns inside 
the shadowed region are not shown. 

scattering length a ID, and one expects to see a similar struc- 
ture. However, one has the advantage of locally tuning the 
interactions using OFR. 

We consider a quasi 2D optical lattice loaded with ^^^Yb 
though other bosonic isotopes will also work. One can in- 
crease U by using OFR and induce a similar structure as 
shown in Fig.[T|b). If one shields a small region of the MI so 
that instead of strong repulsion, bosons in the shielded region 
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interact with each other with a background scattering length 
Qbg = 5.55nm |8|, which should correspond to a small U 
for a lattice with moderate depth. Thus the shielded region 
will be a superfluid. The setup is shown in Fig. [TJd). The 
path labeled (c) on the phase diagram shows how the chem- 
ical potential and interaction strength vary as one goes from 
the shielded region toward the trap edge. The structure of the 
atomic cloud under such an engineered interaction pattern is 
schematically shown in Fig.[TJc). 

It is known that one can generate a superfluid region inside 
a MI core simply by increasing the filling factor. In that way 
the filling factor exceeds 71 = 1 at the trap center so a super- 
fluid core emerges. However, the advantage of using spatially 
varying interactions to induce a superfluid enclave in the MI 
core is that the enclave can, in principle, be induced anywhere 
in the MI core, not just at the trap center We deliberately 
demonstrate this on Fig.[T|c). Moreover, one can induce more 
than one superfluid enclaves in the MI core as long as the res- 
olution allows. With suitably designed masks, one can ex- 
pect that this technique will help realize "printed circuits" of 
atomic devices. We emphasize that similar phenomena can 
be studied using fermions as well. We remark that diffrac- 
tive effects may induce further inhomogeneity effects inside 
the shadowed region such as Arago spot or diffractive oscilla- 
tions. Experimentally one may use a mask with a corrugated 
edge or a combination of masks and lens to minimize diffrac- 
tion. 

So far we consider how to induce a superfluid enclave in- 
side a MI core. One can as well induce a MI enclave inside 
a superfluid. One possible setup is a weakly-interacting su- 
perfluid with a scattering length a^g and filling factor slightly 
above n — 1. Then by imposing OFR in localized regions the 
repulsion U can be increased and when it exceeds the critical 
value, MI enclaves are expected to emerge. Since the den- 
sity profile of a quasi 2D Bose gas in optical lattices can be 
measured quite accurately |20|, one should be able to identify 
those enclaves. In real experiments, however, there should be 
a transient region (or a finite-width wall) between the MI and 
the induced superfluid region. The width is estimated to be of 
order of = 1/ ^/8^mab^, the coherence length of the con- 
densate inside the enclave, because it determines the healing 
length llj. From the parameters shown in Ref. ID, we esti- 
mate that £_c ~ 360nm in the central region, which is reason- 
ably small compare to the Thomas-Fermi radius of that setup, 
which is 5.5/im. 

One future application may be a superfluid-MI-superfluid 
junction constructed by bringing the superfluid enclave close 
to the outer superfluid ring for tunneling phenomena. Such a 
tunneling junction will complement the Josephson junction in 
the absence of optical lattices studied in Ref. Il22ll23l . 



III. INVERTED STRUCTURES OF FERMI GASES 

One major difference between bosonic and fermion systems 
is that in the absence of lattice potentials, the low-temperature 
phase of a Bose gas is always a superfluid. One has to go 
above the BBC temperature to observe a normal phase. To 




Figure 2. (Color online) The trap profiles of n and A (lower inset). 
Solid/dashed lines correspond to spatially varying/uniform interac- 
tions. Here Efq and fc^o are the Fermi energy and momentum at 
the trap center. The proposed experimental setup is shown in the up- 
per inset. The atomic cloud is a unitary Fermi gas induced by OFR, 
but the central region is shielded and remains a normal gas. Possible 
diffraction patterns inside the shadowed region are not shown. 



have coexistence of different phases at low T thus requires 
lattice potentials to induce an insulating (MI) phase. In con- 
trast, for a Fermi gas with attractive interactions, the super- 
fluid transition temperature decreases exponentially with the 
coupling constant in the weakly attractive regime 1 1]. For rea- 
sonably low temperature there is a corresponding critical in- 
teraction strength below which the system is a normal Fermi 
gas. Therefore one can induce interesting structures in har- 
monically trapped Fermi gases without optical lattices. 

One counter-intuitive case of using OFR in Fermi gases is 
to induce a normal-gas region inside a fermionic superfluid. 
In the inset of Figure |2] we show a setup that can generate 
such a structure. It is similar to the setup of Fig.[T](d) except 
there is no optical lattice. We consider Yb with background 
scattering length = — O.lSnm |16| and assume the cloud 
is at unitarity induced by OFR, i.e., a — > —00. Nonetheless, a 
shield blocks OFR at the trap center so the atoms there interact 
with each other via ai,g, which corresponds to extremely small 
attraction ||24l. At low but finite T (say T = O.OSTf, where 
Tp is the local Fermi temperature), the shielded region is a 
normal gas while most of the rest of the cloud is a unitary 
Fermi gas, which is a superfluid. 

For demonstration, we consider the situation where a = 0^ 
for r < Tc and a — —00 for r > rc- We choose Tc = 
{1/4:)Rtf and — Q.lEpQ, where w,. is the radial trap 
frequency and Rtf is the Thomas-Fermi radius. Since the 
unitary limit is around the place where the transition tempera- 
ture shows a maximum in BCS-BEC crossover 1 1 1, this choice 
makes the shell structure easier to be realized in experiments. 
We begin with the ground state close to T = 0. 

In the presence of a harmonic trap, one can use LDA and 
the BCS-Leggett theory to obtain the following equations of 
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states ifTlfTol 

„„.^(,_^), 

ATTh-'a " ^ \2rk ' 2Ek) ■ 

Here Ek = .J (et - fif{r))^ + A{r)^ and = h'^k^/{2m). 
The trap profiles of density and order parameter A of the 
Fermi gas are shown in Figure|2] One immediately sees that in 
order to establish equilibrium with a continuous /ij, the den- 
sity in the superfluid phase has to jump in order to balance 
the diffusion on both sides. Such discontinuities in density 
profiles have been shown in the phase-separated structure of 
polarized Fermi gases due to the same mechanism ifTTl . i.e., 
the balance of /i/, but a normal-gas core cannot be realized if 
the mass of the two species are the same. 

The jump in the density profile serves as direct evidence 
that a superfluid shell encircles a normal-gas core for an un- 
polarized Fermi gas with equal mass. One can also confirm 
the coexistence of a superfluid and a normal gas by measur- 
ing the total energy E of the system from the cloud size f25\. 
The condensation energy from the superfluid at the trap edge 
will make E smaller than that of a cloud consisting of nor- 
mal gas with the same total particle number. To contrast our 
predictions with the case with a uniform interaction, we also 
show in Figure |2]the profiles of n and A of a trapped unitary 
Fermi gas obtained from the same profile of chemical poten- 
tial. There is no jump in the profiles and the whole cloud is a 
superfluid. 

We remark that the sharp discontinuity in the trap profiles of 
n and A is also an artifact of LDA and will be broadened due 
to the kinetic energy terms not considered in LDA. However, 
finite jumps in the density should still be observable based 
on the experimental observations of phase-separated struc- 
tures in polarized Fermi gases 126] . Again we can estimate 
the width of the transient region by the coherence length of 
the superfluid phase, whose magnitude is estimated from the 
BCS formula £,bcs = hvp/A [!}, where vp is the Fermi ve- 
locity. From the parameters of Ref. lfT6l . £,bcs ~ 210nm, 
which is reasonably small compared to the radius of the cloud 
(w 30/im). We emphasize that it is the jump in the density 
that gives an important signature of coexistence of different 
phases and the sharpness of the boundary is less relevant. The 
validity of using LDA and Leggett-BCS theory to describe 
phase-separated structures of polarized Fermi gases 1 11 1 also 
support the validity of using the same method for unpolarized 
Fermi gases with spatially-varying interactions. 

At finite T, the edge of the cloud with low density will be- 
come a normal gas due to thermal excitations. This applies to 
both the case with spatially varying interactions as well as the 
case with a uniform interaction strength. As a consequence, 
a cloud with a uniform interaction has a superfluid core en- 
circled by a normal gas shell at the trap edge. In contrast, 
the structure of a cloud with spatially varying interactions as 
shown in Fig. [2] will be a three-layer one with a normal-gas 
core due to the absence of OFR, a superfluid shell, and an- 
other normal gas shell due to thermal excitations. By intro- 



ducing more complicated patterns of spatially varying interac- 
tions, one may induce structures with more alternating layers, 
which can be identified from jumps in the density profile. 

There is another interesting property of such an inverted 
shell structure. When the cloud is under rotation around the 
center, the superfluid can simply circulate around the normal- 
gas core, which is like a giant vortex core. Since angular mo- 
mentum can be carried in this way, there should be no vortex 
in the superfluid region. The absence of any vortex when the 
cloud is under rotation thus can serve as another evidence that 
the core consists of a normal gas. 

An important contribution of realizing ultra-cold Fermi 
gases with spatially varying interactions is to provide informa- 
tion of how superfluid islands can be stabilized in a quasi-2D 
norma-gas matrix by inhomogeneous interaction strength. For 
fermionic clouds such inhomogeneity can produce modula- 
tions in the density as well as modulations in the order param- 
eter Moreover, the superfluid regions should have higher den- 
sity similar to that shown in Fig.|2]so they are really "islands" 
standing on the normal background. This setup could be use- 
ful in estimating how imbalanced the interactions may be for 
the observed supeconducting islands in underdoped HTSCs as 
shown in Ref. ifTSl . 

When fermions are loaded into optical lattices, more in- 
teresting phases are expected to emerge, including fermionic 
Mott insulator, magnetic phases, rf-wave superfluid, etc |4|. 
Since those phases depend on the strength of interactions, with 
carefully designed masks for OFR laser, one should observe 
many more interesting structures which may not be possible 
in conventional condensed matter systems. 



IV. SUMMARY 

In both bosonic and fermionic systems we present feasi- 
ble experimental setups for generating counter-intuitive struc- 
tures, signatures that can identify the targeted phenomena, 
and experimental techniques that can measure those signa- 
tures. We propose that OFR may enable ultra-cold atomic 
gases to simulate electronics, heavily underdoped cuprate su- 
perconductors, and more. This technique complements oth- 
ers for ultra-cold atoms such as optical lattices, SLM, polar- 
ized Fermi gases, etc. One can envision that combinations 
of those techniques could demonstrate more exciting macro- 
scopic quantum phenomena. 

In this paper we focus on counter-intuitive phenomena in 
equilibrium because they may be realized with minimal ad- 
justment to current experimental setups such as those de- 
scribed in Refs. ||8l [T^ . We point out that spatially varying 
interaction induced by OFR also has great potential in the 
study of quantum transport. For example, one has seen that lo- 
cal chemical potential can be tuned by OFR. A more general 
definition of current is the electrochemical current, which is 
driven by electric field (for charged systems) and/or chemical 
potential. With time-dependent interactions, one can produce 
chemical potential imbalance and a mass current will flow. 
One advantage of using Yb isotopes and OFR for studying 
quantum transport is the tunable interactions between differ- 
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ent isotopes fTS^. A possible direction may be to co-trap dif- through the LANL/LDRD Program, 
ferent species (bosons and fermions) and simulate the sup- 
pression of the current by impurities. 
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